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Abstract 



1 Introduction 



Several proposals have been given as constructive definitions of superstring 
0, 1 1, 1, g Hiii. Type IIB matrix model § |, |, a large 
N reduced model of maximally supersymmetric Yang-Mills theory, is one of 
those proposals. It is defined by the following action: 

S = -L Tr (\[ A ^ Av ] [A ^A v ] + l -i>Y»[A^]\ (1.1) 

where and ip are N x N Hermitian matrices, fi = 1, , , 10 and ip is a 
ten dimensional Majorana-Weyl spinor field. It is formulated in a mani- 
festly covariant way, which is suitable for studying nonperturbative issues of 
superstring theory. Also since it is a simple model of an ensemble of zero di- 
mensional matrices, it is particularly appropriate for numerical simulations. 
In fact we can in principle predict dimensionality of space-time, low energy 
gauge group and matter contents by solving this model. In this paper, we 
pursue this line of analysis and report some results on possible structures 
of the low energy effective theory and the origin of local space-time gauge 
symmetry and diffeomorphism invariance. These analysis gives justification 
of our interpretation of space-time as distributed eigenvalues of the matrices 
A,. 

We first list several important properties of the IIB matrix model. This 
model can be regarded as a large N reduced model of ten dimensional M = 1 
supersymmetric SU(N) Yang-Mills theory. It was shown JT3J that a large 



N gauge theory can be equivalently described by its reduce model, namely 
a model defined on a single point. In this reduction procedure space-time 
translation is represented in the color SU(N) space and eigenvalues of ma- 
trices are interpreted as momenta of fields. As a constructive definition of 
superstring, on the other hand, we will see that we need to interpret eigen- 
values of matrices as coordinates of space-time points, which interpretation 
is T-dual to the above. 

Since our IIB matrix model is defined on a single point, the commutator 
of the supersymmetry 

/ S^A, = zeil> n 9 s 

vanishes up to a field dependent gauge transformation and we can no longer 
interpret this supersymmetry as space-time supersymmetry in the original 
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sense. However, after the reduction, we acquire an extra bosonic symmetry 



SA„ = c M l, (1.3) 



whose transformation is proportional to a unit matrix 1 and an extra super 
symmetry 

h^A, =0 

Combinations of these two supersymmetries (|1.2p and ( |1.4| ) 



1.4) 



qW = qW + g( 2 ), q( 2 ) = i(gw - g( 2 )), (1.5) 

satisfy commutation relations 

[eiQW, e 2 gW] = -2ei 7M e^«J^ (1.6) 



where p^ 1 is the generator of the translation ( |1.3|) and i,j = 1,2. Therefore, 
if we interpret eigenvalues of the matrices A^ as our new space-time coor- 
dinates, the above symmetries can be regarded as ten-dimensional M = 2 
space-time supersymmetries. Since the maximal space-time supersymmetry 
guarantees the existence of graviton, it supports the conjecture that the IIB 
matrix model is a constructive definition of superstring. This is one of the 
major reasons to interpret the eigenvalues of A^ as coordinates of the newly 
emerged space-time. 

The second important and confusing property is that the model has the 



same action as the low energy effective action of D-instantons ||14j| . We 
should emphasize here differences between these two theories since we are 
led to different interpretations of space-time. From the effective theory point 
of view, the eigenvalues represent the coordinates of D-instantons in the 
ten dimensional bulk space-time which we have assumed a priori from the 
beginning of constructing the effective action. On the other hand, from the 
constructive point of view we cannot assume such a bulk space-time itself in 
which matrices live, since not only fields but also the space-time should be 
dynamically generated as a result of dynamics of the matrices and should be 
constructed only from the matrices. The most natural interpretation is that 
the space-time consists of N discretized points and the eigenvalues represent 
their space-time coordinates. One of the main purposes of this paper is to 
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confirm this interpretation of the space-time and make sure that low energy 
excitations propagate in this space-time. 

Final important property is that the type IIB matrix model has no free 
parameters. The coupling constant g can be always absorbed by field redef- 
initions: 

This is reminiscent of string theory where a shift of the string coupling con- 
stant is always absorbed to that of the dilaton vacuum expectation value 
(vev). In the previous analysis of the Schwinger- Dyson equation of IIB ma- 
trix model ||, we have introduced an infrared cut-off e, which gives string 
coupling constant g st = 1/Ne 2 . But through more careful analysis of dy- 
namics of the eigenvalues we have shown that there is no such infrared 
divergences associated with infinitely separated eigenvalues and the infrared 
cutoff e we have introduced by hand can be determined dynamically in terms 
of TV and g. Roughly speaking, as we will study in the present paper, the 
density of the eigenvalues represents dilaton vev and gives the string coupling 
constant. More solid relation requires knowledge of double scaling limit and 
we hope to report in a separate paper. 

There are several reasons we believe that IIB matrix model is a construc- 
tive definition of the type IIB superstring. First this action can be related 
to the Green-Schwarz action of superstring [13 by using the semiclassical 
correspondence in the large N limit: 

- {,}, 

Tr -> / d 2 aJ~g~. (1.8) 



In fact eq. Ql.lD reduces to the Green-Schwarz action in the Schild gauge fllC 



SscMd = / d 2 a[fga{\{X^X v f - % -^{X^ if,}) + pfg\. (1.9) 

Through this correspondence, the eigenvalues of matrices are identified 
with the space-time coordinates X^a) of string world sheet. This is consis- 
tent with our interpretation that the eigenvalues of the matrices represent 
space-time coordinates. 
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The correspondence can go farther beyond the above identification of 
the model with a matrix regularization of the first quantized superstring. 
Namely, we can describe arbitrary number of interacting D-strings and anti- 
D-strings as blocks of matrices, each of which corresponds to matrix regu- 
larization of a string. Off-diagonal blocks induce interactions between these 
strings [I], [L7j . Thus it must be clear that the IIB matrix model is definitely 
not the first quantized theory of a D-string but the full second quantized 
theory. 

It has also been shown that Wilson loops satisfy the string field equa- 
tions of motion for the type IIB superstring in the light cone gauge, which is 
the second evidence for the conjecture that the IIB matrix model is a con- 
structive definition of strings. We consider the following regularized Wilson 

loop |, m 



M 

w{C) = Tr[[] exp{ie{k^ + A n ^)}]. (1.10) 

n=l 

Here k^ are momentum densities distributed along a loop C, and we have 
also introduced fermionic sources A n . e in the argument of the exponential 
is an infrared cutoff associated to space-time extension. In the large N limit, 
e should go to so as to satisfy the double scaling limit. There is still some 
subtly how to take this double scaling limit in which we get an interacting 
string theory. 

Considered as a matrix regularization of the Green Schwartz IIB super- 
string, IIB matrix model describes interacting D-strings. On the other hand 
in the analysis of the Wilson loops IIB matrix model describes joining and 
splitting interactions of fundamental IIB superstrings as Wilson loops. From 
these considerations, it is plausible to conclude that if we can take the correct 
double scaling limit, IIB matrix model becomes a constructive definition of 
the type IIB superstring. Furthermore we believe that all string theories are 
connected by duality transformations, and once we construct a nonperturba- 
tive definition of any one of them, we can describe vacua of any other strings, 
in particular, the true vacuum in which we live. 

Once we are convinced that IIB matrix model is a constructive definition 
of superstring, we can simply go on to do numerical calculations (TPf. But 



before doing so we will clarify in this paper how we can read the space- 
time geometry from matrix configurations and also propose a scenario of 
obtaining local space-time gauge symmetry and diffeomorphism symmetry. 
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It is amusing that these fundamental symmetries can arise from a very simple 
matrix model defined on a single point. 

Dynamics of eigenvalues, that is, dynamical generation of space-time was 
first discussed in our previous paper 0. An effective action of eigenvalues 
can be obtained by integrating all the off-diagonal bosonic and fermionic 
components and then diagonal fermionic coordinates (which we call fermion 
zeromodes). If we quench the bosonic diagonal components x 1 ^ (i = 1...N) 
and neglect the fermion zeromodes £*, the effective action for x^ coincides 
with that of the D = 4 M = 4 supersymmetric Yang Mills theory and van- 
ishes respecting stability of supersymmetric moduli. Inclusions of fermion 
zeromodes and also of non-planar contributions lift the degeneracy and we 
can obtain a nontrivial effective action for space-time dynamics. In our pre- 
vious paper we estimated this effective action by perturbation at one loop, 
which is valid when all eigenvalues are far from one another |x, ; — Xj| >> yjg. 
Of course this one-loop effective action is not sufficient to determine the full 
space-time structure, but we expect that it captures some of the essential 
points for the formation of the space-time. One of important properties of 
the effective action is that as a result of grassmannian integration of the 
fermion zeromodes, space-time points make a network connected locally by 
bond interactions. This becomes important when we extract diffeomorphism 
symmetry from our matrix model. This is discussed in section 4. 

The organization of the paper is as follows. In section 2, we briefly re- 
view our previous results on the dynamics of the space-time and show how 
a network picture of the space-time arises. Here is an analogy with the dy- 
namical triangulation approach to quantum gravity. In section 3 we derive a 
low energy effective theory in some particular eigenvalue distributions con- 
sisting of small clusters of size n. The effective theory acquires local SU(n) 
gauge symmetry and we show that the associated gauge field indeed propa- 
gate in the distribution. This supports our interpretation of the space-time. 
Furthermore we show that a massless fermion field in the adjoint represen- 

3 In this paper we consider matrix configurations near simultaneously diagonalizable 
ones. Then the eigenvalues and the diagonal components of matrices coincide and we 
obtain the same effective action for them at least up to one-loop perturbation in off- 
diagonal components. We expect that the classical space-time picture is valid only in 
such a case. For more general configurations, the classical space-time picture is broken 
and non-commutativity of space-time will become important pcfl . We leave it for future 
analysis. 
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tation appears with a gauge invariant kinetic term. In this way, low energy 
effective theory is formulated as a lattice gauge theory on a dynamically 
generated random lattice. In section 4, we study the origin of gravity and 
diffeomorphism symmetry in the space-time. We show that invariance un- 
der permutations of the eigenvalues leads to the diffeomorphism invariance. 
Background metric is shown to be encoded in the density correlation of the 
eigenvalues. Section 5 is devoted to conclusions and discussions. 



Dynamics of Eigenvalues and Space-time 
Generation 



In this section we briefly review our previous analysis H on the dynamics 
of the eigenvalues to clarify some of important properties. Let us consider 
expansion around the most generic classical moduli where the gauge group 
SU(N) is broken down to [/(l)^ -1 . Then diagonal elements of and if) 
appear as zeromodes while the off-diagonal elements become massive. We 
may hence integrate out the massive modes first and obtain an effective action 
for the diagonal elements. 

We thus decompose A^ into diagonal part X M and off-diagonal part A^. 
We also decompose if) into £ and if>: 



A, 



+ A^; 



x. 



\ 



ip 



v 



(2.1) 



where x l and £^ satisfy the constraints 



N 



1 X fi 



and E^iC* = 0, re- 



spectively, since the gauge group is SU(N). We then integrate out the off- 
diagonal parts A^ and if> and obtain the effective action for supercoordinates 
of space-time S e g[X, £]. The effective action for the space-time coordinates 
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S'cfrfX] can be obtained by further integrating out £: 

J cL4#e _W1 = J dXdie~ s ^ x & 

= J dXe' s ^ x \ (2.2) 

where dX and d^ stand for Ylfs^ II^=o anc ^ Ilili 1 Ila=i dC a i respectively. 
Perturbative expansion in g 2 around diagonal backgrounds (X M ,£) is valid 
when all of the diagonal elements are widely separated from one another: 
\x l — x j \ » g 1 / 2 . At one loop, it is easy to integrate over the off-diagonal 
components. 

After adding a gauge fixing and the Fadeev-Popov ghost term associated 
with the broken symmetry SU(N) /[/(l)* -1 

S s .l + S F . P . = -±Tr([X„A»} 2 ) - ±Tr([X„,b][^,c}), (2.3) 

the action can be expanded up to the second order of the off-diagonal com- 
ponents A^, ip as 

s 2 + s g ,. = ^^((xt-^) 2 ifi^-^V(4-4)^ 

+(C - ovc'.v/ + ~fr»{c - iM.i;/). (2.4) 

The first and the second terms are the kinetic terms for A and ip respectively, 
while the last two terms are A0£ vertices. A bosonic off-diagonal compo- 
nent A l i is transmuted to a fermionic off-diagonal component ip^ emitting 
a fermion zero mode £ l or £ J . This vertex conserves SU(N) indices i and j. 
Note that the propagators for A and ip damp as l/(x l — x j ) 2 or l/(x l — x j ) 
respectively. 

Integration over all the off-diagonal components gives an effective action 
for the zeromodes, x 1 ^ and 

f dAdijdbdc e -(s 2 +s g ,.+s P .P.) = -q det ^ v ^ + s^y 1 

i<j 

= e - s ^°° P ^\ (2.5) 
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where 



C'Pj-. (2.6) 



Here £ y and x^ are abbreviations for — £ J and x^ — x^. 
The effective action can be expanded as 



i<3 

C4 c8 



(ij)' 
C8 

£tr(^ + ^), (2.7) 

which is a sum of all pairs (ij) of space-time points. Here the symbol tr in 
the lower case stands for the trace for Lorentz indices. Other terms in the 
expansion vanish due to the properties of Majorana-Weyl fermions in ten 
dimensions. Note that, since contains two fermion zeromodes, the first 
term Stj\ contains 8 fermion zeromodes £ y along a link (ij) and the second 
term Sf^ contains full 16 fermion zeromodes. Hence, when integrating the 
fermion zeromodes, the second term acts as a delta function for the grassman 
variables <5 16 (£j — Integration over the fermion zeromodes gives the final 
effective action S' e ! fi :loop [X] for the bosonic eigenvalues x*. 

Several comments are in order. The first comment is about supersymme- 
try. The one- loop effective action loop [X, £] has J\f = 2 supersymmetry 
which is a remnant of the original one: 

5x 1 ^ = ieiT^ 1 ,~ q\ 

se =e 2 . 1 • ) 

Transformations for e\ = e 2 and e\ = — e 2 correspond to those generated by 
M = 1 supersymmetry generator Q and its covariant derivative D. In this 
sense zeromodes of x, and & can be viewed as supercoordinates of M = 1 
superspace. 

The next comment is infrared convergence. The partition function is 
obtained by integrating over all the boson and fermion zeromodes. A naive 
expectation is that, since the bosonic coordinates originally parameterize the 
supersymmetric moduli, the integral over x« diverges. However, as proved in 
[H], the integral converges for finite N to all orders of perturbation theory 
due to the damping property of the propagators. At one loop this infrared 
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convergence can be easily seen since two fermion zeromodes in always 
appear with a damping factor l/(x % — x- 7 ) 3 and a condition for saturating 16 
fermion zeromodes requires 1/r 24 damping when two clusters of eigenvalues 
are put apart with distance r. This shows that all points are gathered as a 
single bunch and hence space-time is inseparable. This is consistent with the 
explicit calculations of the partition function [^T| 



The one-loop effective action for X{ is given by further integrating out 
the fermion zeromodes £j. Non-vanishing contributions in the expansion of 
exp(— Slg loop [X, £]) come from terms which saturate all the fermion zero- 
mode integrals d^. Since the terms tr(S^) or tr(S^) in the action S^g 10013 
contains 8 or 16 fermion zeromodes £y respectively along a link (ij), this ex- 
pansion can be visualized as summing over graphs in which some of the atC^ 
pairs among the N space-time points Xi (i = 1...N) are connected by a term 
with 8 fermion zeromodes tr(S^)/A or a term with 16 fermion zeromodes 
(tr(^. ) )) 2 /32 + tr(^. ) )/8: 



J dXd£e- s ^ looPlx >t ] = J dXdi ]T J] 



Gigraph (r/'):bond of G 

„ *•(%)) or AM%\2 , 

4 y 2 y 4 ' 8 

= JdXW[X;G}. (2.9) 

G:graph 

Here H^fXjG] is a Boltzmann weight for a graph G and configuration X. 
In this way, we arrive at network picture of the space-time. The network is 
induced by fermion zeromode integrations. 

Integrations over the fermion zeromodes give interactions along a link (ij) 
of the graph, which is of order l/(x y ) 3 for each fermion zeromode between 
two points (ij). That is, the terms with 8 or 16 fermion zeromodes give bond 
interactions of order l/(x u ) 12 or l/fV- 7 ) 24 ^ In this manner we obtain net- 
works of space-time coordinates connected by bonds with these interactions. 
Since the bond interactions suppress connection between two distant points 



4 Note that the bond interactions induced by terms with 8 fermion zeromodes are not 
scalars and depend on relative angles between several links x l — x J in a graph. On the 
other hand, interactions induced by terms with 16 zeromodes are scalars since acts 
as a delta function for £y . 
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at least by a factor of order l/(x*- J ) 12 , only closer points tend to be connected 
and the network becomes local. Also to saturate the grassman integral, we 
need 16(N — 1) fermion zeromodes. Hence number of bonds is of order N 
which is much smaller than possible number of pairs N(N — 1)/2 in the large 
N limit. In this sense N discretized space-time points are weakly bound. 
All the allowed networks must be summed in determining the dynamics of 
the eigenvalues. This reminds us of summation of all triangulations in the 
dynamical triangulation approach to quantum gravity. We come back to this 
analogy in section 4.. 

As shown in [[|, if we neglect the first type of bond interactions with 
8 fermion zeromodes in (|2.9| ), the £ integration can be performed exactly 
and the dynamics of eigenvalues is governed by the statistical system of a 
branched polymer. Since the Hausdorff dimension of the branched polymer 
is four, we may conjecture that other interactions will make eigenvalue distri- 
bution a smooth four dimensional manifold, which may be classified as a new 
universality class near the branched polymer phase. Numerical simulation is 
still under investigation. 

3 Local Gauge Invariance 

Once we describe the space-time as dynamically generated distribution of the 
eigenvalues, low energy effective theory in the space-time can be obtained by 
solving dynamics of fluctuations around the background X^. Both of the 
space-time and matters are unified in the same matrices and 
should be determined dynamically. Low energy fluctuations are in general 
composites of A^ and ip, and it is natural from the analysis of the Schwinger- 
Dyson equation for the Wilson loops that a local operator in the space-time 
is given by a microscopic limit of the Wilson loop operators, such as 

w (k; O) = Tr[0(A } ip) expiik^A^)}. (3.1) 

Here 0(A,ip) is some operator made of A^ and ip. In order to identify 
the total momentum of this operator as k, the operator 0(A, ip) should be 
invariant under a constant shift of A^, that is, translation in the space-time 
coordinates. 

In the first approximation around the diagonal background X^, the coor- 



10 



dinate representation of this operator is given by 




r d l0 k 



exp(-ik^Xn) w(k;0) 



EO^x-Xi). 



(3.2) 



i=l 



Here we have replaced by X M + and take the leading term. On is the 
ii component of the operator O. Due to the delta function, the operator has 
a support only on the area where eigenvalues distribute. Vanishing of the 
operator w(x; O) outside of the area of the distributed eigenvalues supports 
our interpretation of the space-time. Of course if we take into account higher 
terms in the expansion, the operator becomes dim and extended around the 
background X^. 

We can apply a similar analysis to strings which propagate in the space- 
time. In the 1/N expansion, correlation between Wilson loop operators can 
be evaluated by summing over all surfaces made of Feynman diagrams con- 
necting the Wilson loops at the boundary. This surface is interpreted as 
string world sheet connecting strings at the boundaries. Each SU(N) in- 
dex x 1 ^ of a loop in the diagrams represents a coordinate on the world sheet 
and it takes value in the eigenvalue distribution in the leading approxima- 
tion around the diagonal background X M . Hence string world sheet evolves 
only in the space-time of the eigenvalue distribution and again supports our 
interpretation of the space-time. 

It is generally difficult to obtain how fluctuations propagate in the eigen- 
value distribution, which is reminiscent of the QCD effective theory: Ex- 
citations are expressed as composite operators of microscopic variables and 
their low energy dynamics can be discussed only through the symmetry argu- 
ment, namely the argument based on the chiral symmetry. Also in our case 
we will show that there are eigenvalue distributions around which symme- 
try arguments allow us to discuss low energy dynamics for some excitations. 
Suppose that eigenvalue distribution forms clusters consisting of n eigenval- 
ues. At length scale much larger than the size of each cluster, the SU(N) 
symmetry is broken down to SU(n) m where m = N/n. We can expand A^ 
and ip around such a background X^ similarly to the analysis in the previous 
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section. First write and ip in block forms: 



i, = 



A f 



V 



A 12 
A 22 



?/> 12 
?/> 22 



/l mm 



V 



(3.3) 



Each block AJf or ^ is an n x n matrix and the diagonal blocks can be 
further decomposed 



All 



x 1 ,i + a; 
C 1 + 



(3.4) 



where 1 is an n x n unit matrix and trA" = 0. Here tr means trace for the 
submatrix of n x n. We interpret each cluster of the eigenvalues as a space- 
time point with an internal structure SU(n). Since each SU(n) symmetry 
acts on the variables at the position i independently, the unbroken SU(n) m 
symmetry can be regarded as local gauge symmetry. Indeed under a gauge 
transformation g of the unbroken SU (n) m symmetry 



9 = 



( 9i 



V 



92 



\ 



9m J 



E SU{n) m C SU(N), 



(3.5) 



the diagonal block fields, A™ and ip u , transform as adjoint matters (i.e. site 
variables in the lattice gauge theory) while the off-diagonal block fields, A % i 
and ijfi , as gauge connections (i.e. link variables): 



All 



9i^g\ 



• 9i*P U 9i 

r - 9rt ij gt 



(3.6) 
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Some of the dynamics for low energy excitations is governed by this local 
gauge invariance. Gauge fields live on the links and transform as the link 
variables in the lattice gauge theory. In our case, we have too many such 
fields (at least 10 boson fields A l j for a link (ij) ) but only one unitary link 
variable is assured to be massless by the gauge symmetry and others will 
acquire mass dynamically Therefore, in deriving low energy effective theory, 
we first apply polar decomposition to A l i into unitary and hermitian degrees 
of freedom and identify all the unitary components in various off-diagonal 
block fields by setting them one common field C/ u on each link. 

Among various fields in the dynamically generated random lattice, some 
fields will survive in the low energy effective theory. Here let us consider 
simple fields made of a single microscopic variable or ip. As we will 
see soon, the off-diagonal block fields become massive except the unitary 
component U % K The gauge field U l i is assured to be massless due to the 
gauge invariance. Next, stability of the cluster type eigenvalue distribution 
we have assumed requires that the diagonal bosonic fields A™, which break 
this distribution, should become massive dynamically. Finally, the diagonal 
fermionic block field i\) n can be massless and we should investigate their low 
energy action. 

We henceforth integrate out the off-diagonal blocks first while keeping 
the unitary degrees of freedom and obtain an effective action for other 
variables which survive in the low energy effective theory. 

We add a gauge fixing term for broken generators SU(N)/SU(n) m 

V^Efl^l 2 . (3-7) 

and also its associated Fadeev-Popov term. The action can be expanded as 
before. Lower order terms in the action are composed of the following terms: 

S 2 + S g ,. = ^^tr((x*) 2 \A%\ 2 - fiT^i/,* 

-2$"^ (A^ ji - i/f'A*)). (3.8) 

The first and the second terms are the kinetic terms for the off-diagonal block 
fields A 1 * and ip^ . The third term corresponds to the term of the zeromode 
£ insertion in the previous section. Integrating out the fermion fields in the 
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off-diagonal blocks the third term in ( |3.8| ) becomes 



1 rp^$ 

- £ tr& T^AV-j^iJVjr - tf*J%)). (3.9) 

9 i-Lj \ x 



It consists of a mass-like term proportional to ip T 11 ip u and a hopping term 

ip ip n between i and j. Integration over A 1 ^ can be performed, with its 
unitary component U l i kept fixed, by replacing 

A$ = 



A l Z®Ai l = g*^^U ij ®U ji . (3.10) 



In the second equation, the damping factor l/{x^) 2 corresponds to the prop- 
agation of the hermitian degrees of freedom while the appearance of the link 
variable U 1 ^ corresponds to keeping the unitary degrees of freedom. This 
replacing makes the mass-like term in (|3.9|) vanish and the hopping term 



in 

This is a gauge invariant kinetic term for an adjoint fermion field ?/> M on a 
(random) lattice generated dynamically by the distributed eigenvalues. The 
fields ip u can hop between any pair of points in the space-time, but since the 
hopping parameter is suppressed by l/(x^) 3 , the propagation is expected to 
become local in the continuum limit. 

Other terms in the action generate a plaquette action for the gauge field 
C/ y as follows. Relevant terms in the action are 

S* = - 2 E tr(A^Ai k - AUf)A k jA 1 :. (3.12) 

By integrating out the hermitian degrees of freedom of the off-diagonal blocks 
with the procedure (|3.10Q , this action itself vanishes S4 = 0. However inter- 
actions generated by (64) 2 induce a kinetic term for the gauge field; 

IS? - JC^ (x .. )2( ^. fc) ^ 

(3.13) 
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This is the plaquette action generated by a Wilson loop for the adjoint repre- 
sentation and hence the gauge field C/ u indeed propagates in the space-time 
of eigenvalue distribution. Again the gauge field can hop between any pair of 
space-time points, but the hopping is suppressed by 1/x 8 for distant points 
and we will recover locality in the continuum limit. 

To summarize this section, supposing that distribution of the eigenvalues 
consists of small clusters with size n, we have shown that the low energy effec- 
tive theory contains several massless fields such as the gauge field associated 
with the local SU(n) gauge symmetry and fermion field in the adjoint rep- 
resentation of SU(n) gauge symmetry. Gauge invariant kinetic terms were 
also derived. Thus our system is a lattice gauge theory on a dynamically 
generated random lattice. It is invariant under a permutation for the set of 
the N discrete space-time points, since the permutation group is a sub- 
group of the original SU(N) symmetry. It is the most different point from 
the ordinary lattice gauge theory on a fixed lattice, which becomes impor- 
tant in deriving the diffeomorphism invariance of our model. We will come 
back to this point in the next section. Although the permutation invariance 
requires that all space-time points are equivalent, locality in the space-time 
will be assured due to the suppression of the hopping terms between distant 
points. In general, however, we need a sufficient power for the damping of 
the hopping terms in order to assure locality in the continuum limit. Though 
we do not yet know the real condition for locality, we expect that terms with 
lower powers are canceled due to supersymmetry or by averaging over gauge 



4 Gravity and Diffeomorphism Invariance 



As we have seen in section 2, the one- loop effective action for the space-time 
points is described as a statistical system of iV points whose coordinates are 
x l . Integration over the fermion zeromodes £ gives the Boltzmann weight, 
which depends on a graph (or network) connecting the space-time points 
locally by order N number of bond interactions; 



fields. 




Gigraph 
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W[-X"; G] is a complicated function of a configuration X and a graph G. An 
important property is that the weight is suppressed at least by a damping 
l/(x l — x^) 12 when two points % and j are connected. This system is, of 
course, invariant under permutations Sn Q of N space-time points, which is 
a subgroup of the original symmetry SU(N), while the Boltzmann weight 
for each graph G is not. The invariance is realized by summing over all 
possible graphs. In other words the system becomes permutation invariant by 
rearrangements of bonds in the network of space-time points. This reminds 
us of the dynamical triangulation approach to quantum gravity f22|l , where 



the diffeomorphism invariance is believed to arise from summing all possible 
triangulations. It is amusing that our system satisfies both of locality and 
permutation invariance simultaneously by summing over all possible graphs. 

In this section we see that the permutation invariance of our system ac- 
tually leads to the diffeomorphism invariance. To see how the background 
metric is encoded in the effective action for low energy excitations, let us con- 
sider, as an example, a scalar field <p l propagating in distributed eigenvalues. 
The effective action will be given by 

S = Y: {<pl ~/ )2 f(x* - x?) + £ miff (4.2) 

where f(x) is a function decreasing sufficiently fast at infinity to assure lo- 
cality in the space-time. Introducing the density function of the eigenvalues 

p(x) = Y / 5 {w) (x-x i ) (4.3) 

i 

and a continuous field <p(x) which satisfies (f)(x % ) = l , the action can be 
rewritten as 

S = [ dxdy(p(x)p(y)) ^ ~ f(x -y)+m f dx(p(x))(j)(x) 2 . (4.4) 



Here the expectation (...) for the density and the density correlation means 
that we have taken average over configurations X and networks G of the 



5 In this section we consider general eigenvalue distributions in which all eigenvalues 
have indegenerate space-time coordinates. If we take the cluster type distribution consid- 
ered in the previous section, the permutation symmetry responsible for the diffeomorphism 
invariance should be S^/ n . 
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space-time points. Normalizing the density correlation in terms of the density 



(p(x)p(y)) = (p(x))(p(y))(l + c(x,y)) (4.5) 
and expanding cj)(x) — <f>(y) — (x — y) fl d p '(f)(x) + ■ ■ -, the action becomes 

dy(p(y))(x - y)n{x - y) v f(x - y)(l + c(x,y)) 
d^<p(x)d u <p(x) + m [ dx{p(x))(p(xf ■■■ . (4.6) 



S — — I dx(p(x)) 



2 

2 



This expansion shows that the field 4>(x) propagating in the eigenvalue dis- 
tribution feels the density correlation as the background metric while the 
density itself as vacuum expectation value of the dilaton field. Namely we 
can identify 



9^(x) ~ / dy(p(y))(x-y) (t (x-y) v f(x-y)(l + c(x,y)) (4.7) 
y/ge-*M ~ (p(x)). (4.8) 

If the density correlation respects the original translational and rotational 
symmetry, that is, if they are not spontaneously broken, the metric becomes 
flat g^ u ~ r]^. (Normalization can be absorbed by the dilaton vev.) The fact 
that the background metric is encoded in the density correlations as above 
indicates that our system is general covariant even though the IIB matrix 



model action (|1.1|) defined in flat ten dimensions does not have manifest 
general covariance. 

Then let us see how the diffeomorphism invariance is realized in our 
model. The action (|1 . 1|) is invariant under the permutation Sn of the eigen- 
values, which is a subgroup of SU (N) . Under a permutation 

x i _^ x <i) f or or G s N , (4.9) 
the field 6 l transforms into d> a ^\ Then, from the definition of the continuous 



field 4>(x), we should extend the transformation ( |4.9| ) into x 



x -> ((x) (4.10) 

such that ^(x l ) = x a<y% \ Under this transformation, the eigenvalue density 
transforms clS Qb scalar density and the field 0( scalar field. On the other 
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hand, the metric transforms as a second rank tensor, if the function f(x — y) 
decreases rapidly around x = y and the y integral in (|4.7|) has support only 
near y = x. The tensor property of the metric is also required from the 
invariance of the action under the transformation ( }4.10|) . In this way, the 
invariance under the permutation of the eigenvalues leads to the invariance 
of the low energy effective action under general coordinate transformations. 

Background metric is encoded in the density correlation of the eigen- 
values. Since we have started from the Poincare invariant type IIB matrix 



model action (|1 . 1|) , the density correlation is expected to be translational and 
rotational invariant and we may obtain low energy effective action in a flat 
background. A nontrivial background can be induced dynamically if Lorentz 
symmetry is spontaneously broken and the eigenvalues are non-trivially dis- 
tributed. 

A nontrivial background can be also described by condensing a graviton 
operator Bosonic parts of graviton and dilaton operators are given by 

~ Tr(F, x F\e tkA ) + (fj, <- u) (4.11) 
D(k) ~ Tr(F 2 e ikA ). (4.12) 

Condensation of them induces extra terms in IIB matrix model action; 

5 cond = / rf ME h^(k)S^(k) + h(k)D(k)). (4.13) 

We can similarly obtain an effective action for fluctuations around a diagonal 
background from this modified matrix model action. Condensation of dilaton 
changes the Yang-Mills coupling constant g locally in the space-time. Since 
g is the only dimensionful constant in our model and thus determines the 
fundamental length scale, a local change in g will lead to a local change in 
the eigenvalue density. This is consistent with our earlier discussion that the 
dilaton expectation value is encoded in the eigenvalue density. On the other 
hand, condensation of graviton induces asymmetry of the space-time. For 
condensation of k = graviton mode, it is obvious that the condensation can 
be compensated by a field redefinition of matrices 

^ - (s; + h;)A u , (4.i4) 

and the two models, the original IIB matrix model and the modified one with 
the k = graviton condensation, are directly related through the above field 
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redefinition. The density of the eigenvalues is mapped accordingly and the 
density correlation is expected to become asymmetric in the modified matrix 
model. For more general condensation, if the graviton operator S^ v {x) (co- 
ordinate representation of ( |4.11| )) changes only local property of dynamics of 
the eigenvalues, the density correlation will become asymmetric locally in the 
space-time around x and therefore induces a local change in the background 
metric. 

Our low energy effective action is formulated as a lattice gauge theory 
on a dynamically generated random lattice. Since the lattice itself is gener- 
ated dynamically from matrices, we must sum over all possible graphs. In 
this way, our system is permutation Sn invariant, which is responsible for 
the diffeomorphism invariance. The background metric is encoded in the 
density correlation of the eigenvalues and the low energy effective action be- 
comes manifestly general covariant. The graviton operator is represented 
as fluctuation around the background space-time and constructed from the 
off-diagonal components of the matrices. Microscopic derivation of the prop- 
agation of the graviton is difficult to obtain, but once we have clarified the 
underlying diffeomorphism symmetry, it is natural that the low energy ef- 
fective action for the graviton is described by the Einstein Hilbert action. 
Employing this diffeomorphism invariance and the supersymmetry, we will 
be able to derive the low energy behavior of the graviton multiplets, which 
will be reported in a separate paper. 

5 Conclusions and Discussions 

In this paper we have discussed a possibility to interpret the space-time in 
type IIB matrix model as distributed eigenvalues. We have shown that, if 
we suppose that eigenvalue distribution consists of small clusters with size 
n, the low energy theory acquires SU (N) local space-time gauge symmetry. 
This gauge invariance can predict the existence of the gauge field propagating 
in the space-time of distributed eigenvalues. Also we have obtained a gauge 
invariant kinetic action for a fermion in the adjoint representation of SU(N). 
Low energy behavior for these fields is described in terms of a lattice gauge 
theory on a dynamically generated random lattice and hence supports our 
interpretation of the space-time. 

Since the type IIB matrix model is proposed as a constructive definition 
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of superstring, we need to show the existence of the massless graviton and the 
diffeomorphism invariance of the low energy effective theory. One plausible 
argument for the existence of massless graviton is based on the maximal su- 
persymmetry (D — 10 jV = 2 susy). We have shown that the diffeomorphism 
invariance of our model is originated in the invariance under permutations of 
the eigenvalues. Our model realizes the invariance in an interesting way by 
summing all possible graphs connecting the space-time points. The diffeo- 
morphism invariance restricts the low energy behavior of the model and gives 
another reasoning for the existence of massless graviton. Background metric 
for propagating fields is encoded in the density correlation of the eigenvalues 
while the dilaton vev is encoded in the eigenvalue density. Curved back- 
ground can be described as a nontrivial distribution of eigenvalues whose 
density correlation behaves inhomogeneously. 

Both of these fundamental symmetries, local gauge symmetry and the 
diffeomorphism symmetry, are originated in the SU(N) invariance of the 
matrix model. These symmetries have been attempted to unify by many 
physicists, beginning with Kaluza and Klein [pl| . It is amusing that our 
matrix model can unify them in a natural way and complicated structures 
are hidden in such a simple matrix model. 

There are still several problems unanswered. One is the low energy su- 
persymmetry in the effective theory in a dynamically generated distribution 
of eigenvalues. Since we have obtained a massless gauge boson and an ad- 
joint fermion, we may expect M = 1 low energy supersymmetry. We need to 
clarify how M = 2 supersymmetry in type IIB matrix model acts on the low 
energy fields. Since M = 1 ten dimensional Yang-Mills theory is anomalous 
by itself, other fields (e.g. supergravity multiples) must be required massless 
to cancel anomalies. Also the low energy gauge group might be restricted or 
projected. 

Next problem is locality in the low energy effective theory. The gauge 
invariant kinetic terms ( |3 . 1 1| ) or have hopping parameters suppressed 

due to the inverse powers of distances, but the suppression seems insufficient 
to assure locality or a well behaved infrared property. As for the action 



(|3.11|) , we may obtain more suppressed hopping terms since the link variable 
between distant points can fluctuate much and averaging over the gauge 

field may reduce the strength of effective hopping parameters. 

We have also more fundamental or conceptual issues. We have obtained 

nonabelian gauge symmetry from type IIB matrix model by assuming a par- 
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ticular eigenvalue distribution. This indicates that this vacuum is not a per- 
turbative vacuum of type IIB superstring. Instead we may wonder if this is a 
perturbative vacuum of heterotic string or type I string realized in a nonper- 
turbative way within type IIB matrix model. As we saw in the introduction, 
our matrix model contains both of the world sheets of the fundamental IIB 
string and the D-strings. By a semiclassical correspondence Q1.8Q , we have 
identified IIB superstring in the Schild gauge where tr is interpreted as inte- 
gration over a D-string world sheet. We can also construct an F-string world 
sheet in terms of surfaces made of Feynman diagrams whose SU (N) index 
represents a space-time coordinate of a world sheet point. In both cases, if 
we assume an eigenvalue distribution consisting of small clusters, an internal 
structure appears on the world sheet and hence current algebra may arise. 

Another issue is how to describe global topology in type IIB matrix model. 
The simplest example is a torus compactification. A possible procedure of 
a torus compactification [^J is to identify with + by embedding 
a derivative operator into our matrix configuration. Therefore N is taken 
infinity from the beginning. Since this procedure has a subtlety in the large 
N limit, we need a careful examination of the double scaling limit. 

We also do not yet know how we can describe chiral fermions in lower 
dimensions after compactification. Our compactification procedure looks dif- 
ferent from the ordinary Kaluza-Klein compactification but there is some 
similarity. In the ordinary case, isometry of a compactified space leads to 
local gauge symmetry in four dimensional space-time while in our case the 
internal structure arises from assuming the cluster type eigenvalue distri- 
bution. In our description of four dimensional space-time, eigenvalues are 
confined in four directions and shrunk in the other six directions. Our com- 
pactification therefore can be regarded as a segment compactification in an 
ordinary picture. We need to investigate how chiral fermions can arise in 
such a compactified space-time. 

These problems and issues are under investigations and we want to discuss 
more in the near future. 
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